Finite-size analysis of the Fermi liquid properties of 
the homogeneous electron gas 
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^ ■ Abstract. We analyze the extrapolation to the thermodynamic limit of Fermi liquid properties 

of the homogeneous electron gas in two and three dimensions. Using field theory, we explicitly 
calculate finite-size effects of the total energy, the renormalization factor, and the effective mass 
at the Fermi surface within the random phase approximation (RPA) and discuss the validity 
for general metallic systems. 
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1. Introduction 

We consider the Fermi liquid phase of the (unpolarized) homogeneous electron gas in two and 
three dimensions (d = 2, 3) at zero temperature, T = 0. Quantum Monte Carlo (QMC) 
calculations [H [2] have not only provided the most precise evaluations of the ground state energy, 
E, as a function of density n, but also give access to spectral properties, e.g. the momentum 
distribution, n^, the effective mass, m*, of quasiparticle excitations, and the corresponding 
spectral weight, the renormalization factor Z [21 SI El El El El [TO] . Whereas QMC provides a 
powerful strategy to quantify correlations in Fermi systems, it still suffers from two fundamental 
limitations: the fixed-node error, and the extrapolation to the thermodynamic limit. As shown 
recently, calculations of the Fermi liquid parameters are drastically affected by size effects [H [7] , 
since typical QMC computations are done with N ^ 10 3 electrons, but extrapolation schemes 
based on analytical knowledge of underlying properties significantly reduce the finite size error 
pTJ Q21 [U [7] . Here, to provide further insight, we address the issue of extrapolating calculations 
of finite (small) systems to the infinite volume limit within common language of field theory. 
Strategies of improving the many-body wave function to reduce the fixed-node error in QMC 
are discussed elsewhere [131 HH US] . 

At high densities, r s — > 0, the random phase approximations (RPA) is expected to be valid, 
and we explicitly calculate leading order size corrections within this approximation, and discuss 
their general validity. Here, r s = a/as with a = [2(d — l)Ttn/ d] -1 ^ is the mean inter particle 
distance, as = Ti 2 /me 2 the Bohr radius, and e and m are the charge and mass of the electrons. 
Our results agree with the leading order corrections in [111 01 El [7] where finite size corrections 
have been derived from a formulation directly based on the many-body wave function. In the 
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following, we explicitly concentrate on three fundamental properties, the total energy, E, the 
renormalization factor, Z = n k - — n k +, and the effective mass, m* , at the Fermi surface, kp. 

2. Formalism 

From the dielectric function, e(q,oo), the ground state energy per particle, E, inside a box of 
volume V with periodic boundary conditions, and n = N/V, can be obtained via coupling 
constant integration |16j 



E = E + — E l-v q ~ ~ / d\ / — 
where v q = 2(d — l)-7re 2 /g d_1 is the Coulomb potential and 
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is the ideal gas kinetic energy with E k = h 2 k 2 /2m. In the following we set H = 1. Within the 
RPA [171 US]) the dielectric function is given in terms of the ideal-gas density-density response 
function, xo(k,uj), 

eRPA(k,u) = 1 - v k xo(k,u) (3) 
and the coupling constant integration can be done explicitly 



Erpa = E + — V + / — log e RPj4 (q, iu) 
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Spectral properties are obtained from the single particle Green's function, G(k,u), which 
equals 

G~ l {k,u) = uj + - e k - E(k,uj) (5) 

where the self energy, Ti(k,uj) accounts for interaction effects, and [i is the chemical potential 
which fixes the Fermi energy, Ep = fj,, via Ep = £k F + £(/cF,0;/i = Ep). From a frequency 
integration over the Green's function, the momentum distribution, n&, can be obtained, and the 
magnitude of the jump at the Fermi surface, Z = n, - — n, + , is directly related to the residual 
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of G at (kp,£p). In general, poles of the Green's function give rise to quasi-particle excitations, 
whose energies are given in terms of an effective mass m* at the Fermi surface. 
Within the RPA, the self energy equals [T9ll20l l2T] 

E*PA(fr,w) = -^EjT + ^T^ -Aeiep-e^-d^-s^)] 
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where Ep = k F /2m is the unperturbed Fermi energy. From the self-energy, the renormalization 
factor, Z, at k = kp is given by 



z _ 1 = 1 _as(wo 



and the effective mass, to* , which characterizes the quasi-particle dispersion at the Fermi surface, 
is then obtained via 
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Formally, the RPA contains the leading order corrections beyond the independent particle 
Hartree-Fock approximation, and its validity is restricted to the high density region, r s — > 0. In 
this limit, it is consistent to expand Eq. (|8|). This leads to the so-called on-shell approximation 
(OSA) of the effective mass, 
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In the literature, there exists a long-standing controversy [TH1 E2J [231 IM1 E5], whether Eq. (|BJ) 
or Eq. Q should be used to determine the effective mass. Strictly speaking, the validity of 
the RPA is confined to the high density region where both formulas essentially agree, r s ^ 1, 
and the validity of RPA or RPA-like results at lower densities essentially relies on cancelation 
of errors. However, given the exact self-energy beyond RPA, only Eq. ([8]) provides the effective 
mass, and our finite-size analysis below is based on Eq. ([8]) , since, as we will see, our results are 
not restricted to the RPA. 



3. Thermodynamic limit extrapolation of finite system results 

All formulas of the previous section are valid for systems of finite size, L = V 1 ^, where 
wave vectors are discrete. Perturbative expressions for the Green's function and the dielectric 
functions, as e.g. the RPA expressions, are based on the linked-cluster theorem, valid in a 
grand-canonical ensemble of fixed volume introducing small number fluctuations. Since the non- 
interacting Green's function does not depend explicitly on the system size, the extrapolation from 
the finite to the infinite system corresponds to the conversions of all underlying discrete sums to 
integrals, {2Ti/L) d J2q • • • — > f dq. . ., and finite size corrections are equivalent to discretization 
errors in quadrature. On this basis, the elimination of the size error in the second term on the 
rhs of the total energy expression, Eq. (pQ), is straightforward, leading to the Madelung constant 

^ v q f rfq v q 
vm = Y^ W -J ^y dJ . (10) 

In the following, we show that a similar analysis of the quadrature error can be used to reduce 
the finite size error in other quantities. 

If the integrand is regular, the discretization error is rapidly vanishing (of order 1/N or 
higher). The main contribution of order N~ a with a < 1 comes from non-analytical points in 
the integrand. From the general structure of the perturbation expansion there are only certain 
discrete points where non-analyticities can develop: in the long wavelength limit, k = 0, where 
the Coulomb potential is singular, and for integer multiples of the Fermi wave- vector k = j kp 
(j = 1,2,...) due to the sharp Fermi surface of the underlying wave function. 

Explicitly, the sharp Fermi surface at kp leads to well-known shell oscillations in the energy, 
already present in the energy of the ideal Fermi gas, Eq. ([2]). Twisted boundary conditions 
strongly reduce these oscillations [26], and grand-canonical twist averaging (GTABC) completely 
eliminates them [TTJ. Within GTABC, we have 

e gtabc = — f — V" 6{e F - e k+ ^)e k+ ^i = - [ — , d 0{e F - e k )e k (11) 



and Eq TABC = E (N ->• oo). Similarly, GTABC also eliminates size effects of other non- 
interacting properties, e.g. in the ideal gas density response, Xo TABC {q, w ) = Xo(q,w; N — > oo), 
as can be seen from the spectral representation 

aTABC(„ , ,\ _ 2 f d $ Q( £ F ~ gk+g) ~ 8{£f ~ gk+g+q) no v 

V J|0 Q |<7r/L (2vr/L) d ka= ^ m /L w + £k+,? ~ £ k+tf+ q + m? 

although xo is still only given at discrete wave vectors (q a = 2-irj/L with integer j). Therefore, 
twisted boundary conditions as well as GTABC do not modify the discrete sum in the exchange- 
correlation energy on the rhs of Eq. (pQ). However, GATBC lead to an accelerated convergence 
of the integrand towards the thermodynamic limit on all allowed wave vectors. 

Within GTABC, we assume that size effects in the integrands of all expressions can be 
neglected, and focus on the presumably largest contributions from the Coulomb singularity at 
k = 0. For the total energy, Eq. ([T]), we expect that the leading order size corrections of the 
exchange-correlation energy, 5E, is given by 
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Similarly, the effect of the thermodynamic limit extrapolation on the renormalization factor and 
the effective mass can be obtained from explicit (approximate) expressions of the self-energy, 
Eq. ([U]), together with Eq. (0) and Eq. ©. 

4. Results 

The Lindhard function, xoiQi^)-, has the following limiting form in the long wave length limit 

Xo(k,u) re nk 2 /muj 2 , k ->■ (14) 

and we have 

e(k,u)) « 1 -uj 2 p {k)/uo 2 } k^O (15) 
where the plasma frequency is given by u p {k) = {nvuk 2 /to) 1 / 2 . 

4-1. Total energy 

Inserting the limiting form of the dielectric function in Eq. (|H), we obtain 

n J\<ia\<Tr/L (2TT) a J 27T y 11 J|q Q \< W /L {2lT) a 2 

and size corrections in the energy have a simple interpretation in terms of zero point energy of 
long wavelength plasmon modes. Explicitly, in three dimensions, we have 

h 2 /~3~ 1 

5Em= 2^J^N {l7) 

in agreement with Ref. |llj . and 



«W = ,f(4r)V* 1 ^ (18) 
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in two dimensions [I] where S2 = £^ iig a |<i \ a \ 1 ^ 2 = n Jo^ n (l + u 2 ) 1 ^ 4 ~ 1.362073657 accounts 
for the cubic geometry. 



4-2. Renormalization factor and effective mass 

For the calculation of the renormalization factor and the effective mass, we need the derivatives 
of the self energy at the Fermi surface. Within the RPA, we have 
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Since Xo(<?> 0) = —dn/e F for q — > 0, the static dielectric function diverges in the long wave length 
limit, 



CRPA(q, 0) = 1 + dnv q /e F ~ v q ~ q 
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providing perfect screening, and the dominant size correction, A, of dTi/duj exactly equals that 
of —dT,/dk, and is given by 
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in three dimension [7] with c 3 = ^ jj K1 |g|- 2 = £ jjdu ln [^ 2) ~ 1.221374804, and 
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-3/2 _ i /^(l + u 2 )- 3 / A ~ 1.057929920. The 
leading order size corrections for the renormalization factor are then given by 



in two dimensions [3] with C2 = £^ Ju a |<i 



and, using Eq. ([8]), the effective mass correction writes 
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Note, that, within the validity of the RPA, we may put Z = 1 on the rhs of Eq. (|28|) . The 
slow decay of A with increasing system size can lead to important qualitative and quantitative 
changes of spectral quantities when properly extrapolated to the thermodynamic limit with 
respect to calculations which rely on linear extrapolation schemes. 



5. Discussion 

We have discussed size effects of the total energy, the renormalization factor, and the effective 
mass. Particularly, in two dimensions, Z and m* show an unexpected slowly decaying size 
effects ~ iV -1 / 4 . In general, we expect that these size effects are also present in other response 
quantities, e.g. the spin susceptibility, as they are connected to m* by relations based on 
Landau's Fermi liquid theory. Our results on the finite size corrections are quite general and 
do not rely on the validity of the RPA, as they are based on the large wave length behavior of 
the dielectric function, Eq. (|15|) . which applies to all metallic systems. Further, all results agree 
with the leading order corrections in [H [7] where finite size corrections have been derived 
from a formulation directly based on the many-body wave function. 
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